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We introduce a general framework for analyzing the thermodynamics of small systems that are 
driven by both a periodic temperature variation and some external parameter modulating their 
energy. This set-up covers, in particular, periodic micro and nano-heat engines. In a first step, we 
show how to express total entropy production by properly identified time-independent affinities and 
currents without making a linear response assumption. In linear response, kinetic coefficients akin 
to Onsager coefficients can be identified. Specializing to a Fokker-Planck type dynamics, we show 
that these coefficients can be expressed as a sum of an adiabatic contribution and one reminiscent of 
a Green-Kubo expression that contains deviations from adiabaticity. Furthermore, we show that the 
generalized kinetic coefficients fulfill an Onsager-Casimir type symmetry tracing back to microscopic 
reversibility. This symmetry allows for non-identical off-diagonal coefficients if the driving protocols 
are not symmetric under time-reversal. We then derive a novel constraint on the kinetic coefficients 
that is sharper than the second law and provides an efficiency-dependent bound on power. As 
one consequence, we can prove that the power vanishes at least linearly when approaching Carnot 
efficiency. We illustrate our general framework by explicitly working out the paradigmatic case of 
a Brownian heat engine realized by a colloidal particle in a time-dependent harmonic trap subject 
to a periodic temperature profile. This case study reveals inter alia that our new general bound on 
power is asymptotically tight. 

PACS numbers: 05.40.-a, 05.70.Ln 


I. INTRODUCTION 

Thermodynamic processes on the micro- and nano¬ 
scale in systems driven out of equilibrium by periodically 
changing control parameters like an external force or the 
temperature of their environment, can be scrutinized un¬ 
der the microscope by virtue of precise measurements 
of characteristic quantities such as applied work or ex¬ 
changed heat [ll-Q. Despite such groundbreaking results, 
so far, no general theoretical framework for the thermo¬ 
dynamic description of periodically driven systems be¬ 
yond the quasi-static limit is available. 

In principle, the concepts of irreversible thermody¬ 
namics, a phenomenological but powerful theory, which, 
building on the principle of local equilibrium, furnishes 
non-equilibrium stea^ states with a universal thermo¬ 
dynamics structure [7|, can be transferred to periodic 
systems These results are, however, crucially tied 

to specific models and require rather involved and non- 
intuitive definitions for currents and affinities. 

In this paper, we overcome these limitations. Start¬ 
ing from the first law formulated for an arbitrary system 
driven out of equilibrium by both, periodic perturbations 
of its Hamiltonian and the temperature of a surrounding 
heat bath, we obtain natural and perfectly general iden¬ 
tifications of fluxes and affinities. Since they are defined 
on the level of cycle averages, these quantities are time- 
independent. Nevertheless, our new formalism captures 
essential finite-time properties of the driven system and 
permits a discussion of quantities like power, which are 
out of reach for the laws of classical thermodynamics due 


to its notorious lack of time scales. 

Moreover, our approach provides a universal prescrip¬ 
tion for the construction of kinetic coefficients, which 
fully characterize the system in the linear response 
regime. By using a well established and rather general 
stochastic approach based on a Fokker-Planck equation 
to describe the underlying time-dependent dynamics [HI, 
we prove that these coefficients are interrelated by a re¬ 
markable symmetry, whi ch, j ust like Onsager’s celebrated 
reciprocity relations [1^ [l^ , can be traced back to mi¬ 
croscopic reversibility. 

Ever since James Watt’s steam engine, the urge to ex¬ 
plore the fundamental principles governing the perfor¬ 
mance of cyclic heat engines, was one of the major quests 
in thermodynamics. Two key figures of merit are partic¬ 
ularly important in this context. While efficiency, the 
first one, is universally bounded by the Carnot value as a 
direct consequence of the second law, similar constraints 
on power, the second one, could, so far, be obtained only 
within specific setups, see for example [l^ - [l8l| . 

Thermoelectric heat engines, which, in contrast to 
cyclic ones, work in a steady state [T^, [2^ are like¬ 
wise subject of substantial research efforts concerning ef¬ 
ficiency [2l|, , power and the relation between 

these quantities [2^, [2^ . At least in the linear response 
regime, however, much more general results are currently 
available for this type of heat engines than for their peri¬ 
odic counterparts. Specifically, it can be shown that, in 
the linear regime, the power of such devices is bounded 
by a quadratic function of efficiency that vanishes at the 
Carnot value rjc and attains its maximum at the Curzon- 
Ahlborn value ??c/2 This result follows from a 
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quite general analysis within the framework of linear ir¬ 
reversible thermodynamics. It must, however, be recon¬ 
sidered in the presence of a magnetic field, which breaks 
the Onsager symmetry, an issue which is currently under 
active discussion jlO) f3ll - l40l| . 

Here, by applying our new formalism, we establish the 
aforementioned quadratic bound on power for cyclic heat 
engines in linear response. In particular, by using a novel 
method, which does not require any additional assump¬ 
tions, we prove that this bound still holds if the matrix of 
kinetic coefficients is not symmetric, which is the generic 
case for periodically driven systems. We emphasize that, 
by now, an analogous result for thermoelectric heat en¬ 
gines is not available on such a general level. To complete 
our analysis, we show that this bound on power is tight 
within a paradigmatic model of a Brownian heat engine, 
which was originally proposed in 0 and recently real¬ 
ized in a landmark experiment Q. 

The rest of this paper is structured in three major 
parts. In section II, we develop our general formalism 
and prove a generalized reciprocity relation for the kinetic 
coefficients of periodic systems. Section III is devoted to 
the discussion of cyclic heat engines in the linear response 
regime and the derivation of a general bound on power. 
We illustrate our results by considering a simple model 
system in section IV. Finally, we conclude in section V. 

II. FRAMEWORK 

In this section, we will demonstrate that the notions 
of irreversible thermodynamics can be transferred to pe¬ 
riodically driven systems. 

A. Nonlinear Regime 

We begin with a brief review of the energetics of driven 
systems in thermal contact with a heat bath [s^ . Specif¬ 
ically, we consider a classical system with degrees of free¬ 
dom X E (xi,... ,Xn) and time-dependent Hamiltonian 

i7(x,t) E i7o(x) H-AiJg™(x,t), (1) 



X 


FIG. 1. Non-equilibrium periodic state of a system with one 
degree of freedom x in a sinusoidally shifted harmonic poten¬ 
tial, represented by gray parabolas, which is embedded in an 
environment with periodically changing temperature as indi¬ 
cated by the periodic color gradient. The solid line in the x-t 
plane shows the motion of the center of the Gaussian phase 
space distribution, which lags behind position of the potential 
minimum shown as dashed line. At any time t, the width of 
the colored region equals the width of the phase space distri¬ 
bution, which varies according to the temperature. 

the system picks up the heat 

Qit) = JcT'x (4) 

from the environment as stipulated by the first law 

U{t)^Q{t)-W{t). (5) 

We will now pass from time-dependent to constant 
variables by exploiting the periodic boundary conditions 

iJ(x,t-tr) = i?(x,t) and T(t + T) =T(t), (6) 


which is immersed in a heat bath, whose temperature 
T{t) oscillates between the two values and > 
Tc- Here, denotes an externally controlled di¬ 

mensionless function of order 1 and AH quantifies the 
strength of this time-dependent perturbation. The power 
extracted by the controller thus reads 

W{t) = -j(ryiH{yi,t)p{yL,t), (2) 

where p(x, t) denotes the normalized probability density 
to find the system in the state x at the time t and dots 
indicate time derivatives throughout the paper. To com¬ 
pensate for this loss in internal energy 

U{t) = y"d"x i7(x, t)p(x, t), (3) 


where T is the length of one operation cycle. Quite nat¬ 
urally, we invoke the assumption that, given these condi¬ 
tions, the time evolution of the probability density p(x, t) 
eventually converges to a periodic limit 


p''{x,t) = p''{x,t + T) (7) 

as illustrated in Fig. [T] for a simple system. Once this 
periodic state is reached, the average entropy production 
per cycle arises solely due to heat exchange between the 
system and the environment, since, due to the period¬ 
icity of the the distribution p‘^(x,t), no net entropy is 
generated in the system in a full cycle, i.e., we have 


S=- 


1 

r 



Qjt) 

Tit)- 


( 8 ) 
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By inserting (jH) into dS]) and parameterizing T{t) by a 
dimensionless function 0 < 7 q(t) < 1 such that 

T(t) = _ TcT], _ 

^rh + (r,-Th) 7 ,W’ 

it is straightforward to derive the expression 


7~ T 1 ^ ' j ^ 

S=—— dtj d"xg^(x,t)/(x,t) 

+ ( 10 ) 

using one integration by parts with respect to time. The 

corresponding boundary terms vanish due to the period¬ 
icity of the involved quantities. Obviously (HOD can be 
cast in the generic form Q 

S = TqJq, (11) 

by identifying the work flux 

Jw = Y J^dtJ(rxg^{x,t)p''{x,t), (12) 

the generalized heat flux 

Jq=^j^dt j d^x-fq{t)H{x,t)if’{x,t) (13) 

and the affinities 


T^ = AHIT, and .F, e 1/7; - l/Th. (14) 

Within only a few lines, we have thus obtained our first 
main result, namely, we recovered for periodically time- 
dependent systems the structure of irreversible thermo¬ 
dynamics. The key point here is the identification of ap¬ 
propriate pairs of affinities and fluxes, whose respective 
products sum up to the total entropy production. 

For later purposes, we note that, after one integration 
by part with respect to t, the heat flux m can be rewrit¬ 
ten as 


Jg = ^ dty"d”x gqp'^{x,t) 

T T ^ 'I ^ 

+ — j^dtj d'^x-iq{t)g^{x,t)f{x,t) (15) 

where 

5 g(x,t) E-i7o(x)7,(t). (16) 


B. Linear Response Regime 

1. Kinetic Coefficients 


with linearized affinities 

ffi^ = AHlT^ and Tq = ATjT^ + (D {A^) (18) 

and kinetic coefficients 
B J 

Lc/ 3 =^:^ for a,j3 = w,q. (19) 

The entropy production ([5D thus reduces to 

5= ^ L^pVp. (20) 

06,/3=w,g 

To guarantee that this expression is nonnegative for any 
Ta as stipulated by the second law, the kinetic coeffi¬ 
cients must obey the constraints 

d^wwi d^qq — and Lyj'^jjLqq — (^Lyjq + Lqyj )V4>0, (21) 

which we prove explicitly in Sec. 1111 BI for a large class of 
systems. It is, however, not evident at this stage whether 
a reciprocity relation relating L^p with Lp^ or any fur¬ 
ther constraints exist. 


2. Adiabatic Limit 

As a first step, we investigate the adiabatic regime, 
which is characterized by the Hamiltonian H(x,t) and 
the temperature T(t) changing slowly enough in time 
such that the system effectively passes through a se¬ 
quence of equilibrium states, i.e., 

/(x,t) =exp[-H{x,t)/(kBT(t))]/Z{t) (22) 


with 


Z{t) E Jcffix exp[-H{x,t)/{kBT{t))] (23) 

and kB denoting Boltmann’s constant. Expanding (12^ 
to linear order in AH and AT and inserting the result 
into (USD, (Ha and dUD gives the universal expression 

Ll^p = ~{{S9a6gp)) (24) 

for the adiabatic kinetic coefficients. Here, we introduced 
the notations 

((A)) E 1 Jjdt {A{t)) E 1 Jjdt fd-xA{x,t)p^^x) 

(25) 

and 


We now focus on the linear regime with respect to the 
temporal gradients AH and AT = T^^-T^. By expanding 
the fluxes (TT^ and (fT^ . we obtain 


dw ~ I^WW'd' w LwqdFq + 0{A^'), 
Jq = Lqwdf^w + LqqdFq H- 0{A^) 


( 17 ) 


dA(x, t) = A(x, t) - (A(x, t)) (26) 

for any quantity A(x, t) and the equilibrium distribution 
of the unperturbed system 

p^'^Cx) E exp [-77o(x)/(fcsrc)] /Zq 


(27) 
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with Zq denoting the canonical partition function. 

Notably, the coefficients (IM)) are fully antisymmetric, 
i.e., 

rad _ rad /'OS'l 

^al3 - ~^l3a- 

As might be expected, this property, which can be proven 
by a simple integration by parts with respect to f, im¬ 
plies vanishing entropy production (I20|) in the adiabatic 
limit. This avoidance of dissipation can, however, be only 
achieved for an infinite cycle duration T and therefore in¬ 
evitably comes with vanishing fluxes Ja- 


The dagger showing up in (l32l) . from here onwards, des¬ 
ignates the adjoint of the respective operator. Note that, 
while in (1551) the operator L°(x) acts on the function 
p®‘^(x), (1551) is to be read as an operator identity becom¬ 
ing meaningful when applied to a specific function of x. 
Physically, the relation (1551) means that, once the system 
has reached its equilibrium state, the rate of transitions 
from the state x to the state x' is balanced by the rate 
of transitions in the reverse direction. 

The equilibrium Fokker-Planck operator L*^(x) can be 
naturally decomposed in a reversible and an irreversible 
contribution 


3. Stochastic Dynamics 


For further investigations of the kinetic coefficients, we 
have to specify the dynamics, which governs the time 
evolution of the probability density p(x, <). Having in 
mind, in particular, mesoscopic systems surrounded by 
a fluctuating medium, a suitable choice is given by the 
Fokker-Planck equation [4l| 

dtp{yi, t) = L(x, t)p{x, t) (29) 


with 


L(x,t) = -d^,D,ix,H,T) + (30) 


where summation over identical indices is understood and 
natural boundary conditions are assumed. The Hamilto¬ 
nian H (x, t) and the temperature T(t) enter via the drift 
and diffusion coefficients Di{x,H,T) and Dnix H,T), 
which thus become implicitly time-dependent 42|. 

We will now formulate a set of conditions on the gen¬ 
eral Fokker-Planck operator (15(11) to adapt it to the phys¬ 
ical situation that we wish to discuss here. Since micro- 
reversibility plays a crucial role in linear irreversible ther¬ 
modynamics, we have to ensure that our theory complies 
with this fundamental principle. To this end, first, at any 
time t any possible state x of the system must be asso¬ 
ciated with the same energy as the time-reversed state 
£x E {siXi,... ,£nXn) with £i = 1 for even and Si = -1 for 
odd variables, i.e.. 


H{x,t) = (31) 

where, throughout the paper, the transformation x ^ ex 
is meant to include the reversal of external magnetic 
fields. Second, the unperturbed Fokker-Planck opera¬ 
tor L°(x) E L(x,t)|^^Q must obey the detailed balance 
condition (4ll| 


L°(x)p'^'5(x)=p'='i(x)L°t(ex), (32) 

for the canonical distribution (I27|) . which uniquely satis¬ 
fies 


(33) 


LL(x)E(L°(x)-L°(ex))/2, (34) 

L°„(x)E(L°(x) + L°(ex))/2, (35) 

which are characterized by their respective behavior un¬ 
der time reversal. While the irreversible part accounts 
for dissipative effects induced by the presence of the heat 
bath, the reversible part describes the intrinsic coupling 
of the system’s degrees of freedom, which is not directly 
affected by the fluctuating environment. Since this au¬ 
tonomous part of the dynamics should preserve the inter¬ 
nal energy of the system, we have to impose the condition 

l°Ux)Ho{x) = 0. (36) 

We note that this consideration does not play a role in the 
overdamped limit, within which the entire time evolution 
of the system is effectively irreversible due to strongly 
dominating friction forces. 

The notion of detailed balance can not be immediately 
generalized to situations with external driving and time- 
dependent temperature. However, in analogy with (1331) . 
the full Fokker-Planck operator L(x, t) can still be char¬ 
acterized by the weaker property 

L(x,f)exp[-7J(x,f)/(fcB'r(i))] = 0, (37) 

which is naturally obeyed in the absence of nonconser¬ 
vative forces and guarantees that the system follows the 
correct thermal equilibrium state if the Hamiltonian and 
the temperature are varied quasi-statically. 


4- Finite-Time Kinetic Coefficients 

In the linear response regime, the Fokker-Planck op¬ 
erator L(t) showing up in (1551) can be replaced by the 
expansion 

L(f) E L° -H H- ATL^(t) s- (38) 

where, for simplicity, from (1381) onwards, we notationally 
suppress the dependence of any operator on x, whenever 
there is no need to indicate it explicitly. The Fokker- 
Planck equation (1551) can then be solved with due consid¬ 
eration of the boundary condition dS]) by treating L^{t) 


L°(x)p"'’(x) = 0. 
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and L^(i) as first order perturbations. The result of this to the compact expression 
standard procedure [4l|, reads 


r' oo Q 

+ X; AX / dre*-a^(t-r)p®'i(x)+0(A2). (39) 

X=H,T -^0 

After some algebra involving condition (I37|) , which we rel¬ 
egate to appendix [A1 for convenience, this solution leads 


1 

Lap = L^ap + ^ Jo ((<55a(0);%(-T))) (40) 


for the kinetic coefficients, where the generalized equilib¬ 
rium correlation function is defined as 



A(x, <1 -H t)e'-“(‘i-*=)B(x, t 2 + <)p®'^(x) 
B(x, t 2 + t)e'-“(*^-*Azl(x, ti + t)p^^{x) 


for (<i > t 2 ) 
for {ti < t 2 ) 


(41) 


for arbitrary quantities A(x, t) and We recall 

the definition (1^51) of the d-notation. 

The expression (1401) admits an illuminating physical 
interpretation. It shows that the kinetic coefficients of 
periodically driven can be decomposed into an adiabatic 
contribution independently identified in (1241) and a finite¬ 
time correction, which has the form of an equilibrium 
correlation function. This result might therefore be re¬ 
garded as a generalization of the well established Green- 
Kubo relations, which relate linear transport coefficients 
like electric or thermal conductivity to equilibrium cor¬ 
relation functions of the corresponding currents [d^. In 
our case, the role of the currents is played by the fluc¬ 
tuation variables <5 Jq(x, t) and the ensemble average is 
augmented by a temporal average over one operation cy¬ 
cle. We note that a similar expression has been obtained 
for the special case of the effective diffusion constant of 
a periodically rocked Brownian motor in . 


5. Reciprocity Relations 

Time-reversal symmetry of microscopic dynamics ap¬ 
pears as the detailed balance condition on the level of the 
Fokker-Planck equation (1321) . By using this relation and 
the Green-Kubo type formula (|4ni) . it is straightforward 
to derive the generalized reciprocity relation 

Lap = Lpa , 

(42) 

where, the Onsager coefficients are considered as func¬ 
tionals of the time-dependent Hamiltonian and tempera¬ 
ture and an external magnetic field B. The technical de¬ 
tails of the derivation leading to the relation (l42l) , which 
constitutes our second main result, can be found in ap¬ 
pendix m Here, we emphasize that, although the setup 
of this paper differs significantly from the one Onsager 
dealt with in his pioneering work the symmetry 

(HU), which constitutes our second main result, and the 
original Onsager relations share microscopic reversibility 
as the common physical origin. Since, in the presence of 


time-dependent driving, full time reversal especially in¬ 
cludes reversal of the driving protocols, naturally, these 
reversed protocols show up in (1^ . 

The symmetry (1421) holds individually for both, the 
adiabatic kinetic coefficients and the finite-time correc¬ 
tion showing up in (HOI) . Given the general relation (1^ , 
it follows that the must vanish if the driving proto¬ 
cols are symmetric under time-reversal. 

The additional relation 

Lap['yw{t),^q{t),B] = Lpa [7g(t),7u,(t),-B], (43) 

can be proven if the driving introduced in ([T]) 

factorizes according to 

= gw{x)^w{t). (44) 

Hence, the off-diagonal kinetic coefficients change place if 
the magnetic field is reversed and the respective protocols 
determining the time dependence of the Hamiltonian and 
the temperature are interchanged. This symmetry does 
not involve the reversed protocols. It is, however, less 
universal than (1^ , since it requires the special structure 
(|3H), see appendix IB] for details. 


III. CYCLIC STOCHASTIC HEAT ENGINES IN 
LINEAR RESPONSE 

As a key application of our new approach we will dis¬ 
cuss the performance of stochastic heat engines. 


A. Power and Efficiency 


The main two benchmark parameters here are, first, 
the average power output per operation cycle 

P=-^ dt J d”x H(x,f)p‘'(x,t) =-TcPwJw (45) 

and, second, the efficiency 

g = PlJq = -TcPwJwlJq, 


(46) 
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which is bounded by the Carnot value rye = 1 “ as 

a direct consequence of the second law 5 > 0. The latter 
figure, which is naturally suggested by the representa¬ 
tion CD of the entropy production per cycle, should be 
regarded as a generalization of the conventional thermo¬ 
dynamic efficiency defined for a heat engine operating 
between two reservoirs of respectively constant tempera¬ 
ture. Our formalism includes this scenario as the special 
case where 7 q(t) is chosen as a step function 


7?(i) = 


for 0 < t < 7i 
for Ti<t <T 


(47) 


with 0 < 7i < T such that the system is in contact with 
the hot temperature Th in the first part of the cycle and 
the cold temperature Tc in the second one. 

Under linear response conditions, which we will assume 
for the rest of this section, the fluxes Ja can be eliminated 
using (HZl such that the expressions (H51) and (H51) reduce 
to 


P = -TePw {LwwJ^w + Lwqpq) 


(48) 


and 


77 = 


TcPwjLww^w + LuiqTq) 

Lqw^w LqqPq 


(49) 


respectively. Clearly, these figures are crucially deter¬ 
mined by the kinetic coefficients La/s- In contrast to 
the thermoelectric case, where the reciprocity relation 
La /3 = LPa holds without magnetic fields, the analysis of 
the preceding section has revealed that for cyclic heat en¬ 
gines this symmetry is typically broken if the driving pro¬ 
tocols are not invariant under time reversal. As pointed 
out by Benenti et al. , a non symmetric matrix of ki¬ 
netic coefficients leads to profound consequences for the 
performance of thermoelectric devices including the op¬ 
tion of Carnot efficiency at finite power. 

These results apply similarly to the systems con¬ 
sidered here, since our theoretical framework is struc¬ 
turally equivalent to the standard theory of linear irre¬ 
versible thermodynamics used as a starting point in [sij . 
To demonstrate this correspondence explicitly, following 
31|, we define the dimensionless parameters 


_ Pwq 
Pqw 


and 


y = 


PwqPqw 


Pww Pqq 


T T 

^wq^qw 


(50) 


which, due to the second law (EH), are related by the 
inequalities 


h{x) <j/<0fora;<0, 0<y< h{x) for cc > 0 (51) 


with 

{x - 1)^ 

By optimizing (l4^ and (j49]), respectively, as functions of 
Pw the expressions 

?7max(a;,?/) = (53) 

V y + 1 + 1 



VMc 

FIG. 2. Plots of the maximum power (Ezl) in units of Po = 
TcJ-qLqq as a function of the normalized efficiency fj = rjjric 
for selected values of the asymmetry parameter a; > 1 in the 
upper and a; < -1 in the lower panel. For x ->■ ±00 , convergence 
to PjPo = fj is observed. The dotted black lines correspond to 
the maxima of dSH for 0 < a; < 00 in the first and -00 < a; < -1 
in the second plot thus indicating the relation between maxi¬ 
mum power and efficiency at maximum power. The apparent 
divergence for negative x occurs in the limit a; -> -1 [^ . 


for maximum efficiency and 


V*{x,y) (54) 

4+ 2y 


for efficiency at maximum power [30l| are obtained, where 
rjc ~ at /Tq = TePq denotes the Carnot efficiency in the 
linear regime. For y = h(x) and |a:| > 0, the maximum 
efficiency equals rjc and the efficiency at maximum power 
can exceed the Curzon-Ahlborn value 7?c/2 IdSl . 

1 ^ . reaching even rjc in the limit a; ^ ± 00 . 

For a quantitative assessment of the relation between 
power and efficiency, following , we consider the max- 
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imum power at given efficiency 


P{V,X,y) = TcTgLqqtJ 


( x(2 + y)-yr] 
\ 2x{l + y) 


+ 


y‘^(x + r])^ yr] 
4x^(1+ 2 /)^ x{l + y) 


(55) 


right 2x2- submatrix, we recover the inequality (I21L 
which we inferred from the second law on the phenomeno¬ 
logical level and now have proven explicitly. Second, by 
evaluating the determinant of A, we get the new con¬ 
straint 


Lqq < Nqq 


PwwPqq 


(A 


wq 


+ Lqyj)^ jA 


LwwLqq PwqPqw 


(62) 


as a joint benchmark parameter, which is found by elim¬ 
inating Au, in (1451) in favor of 


T] = r]/r]c (56) 

using dm). To ensure that the Carnot value rjc is in¬ 
cluded in the range of accessible efficiencies, y must be 
replaced by its bound h{x) in (1551) . The resulting func¬ 
tion 


P{r],x) = P{r],x,h{x)) (57) 

is plotted in Fig. O While P{t],x) reduces to 

P{rj,l)=nP^^Lqqfjil-fj) (58) 

in the symmetric case x = 1 and thus vanishes linearly 
in the limit ry ->• rye, strikingly, we observe that Carnot 
efficiency might be reached at finite power (4^ 

P{r]c,x) = TcPgLqq- — \ for |x| > 1. (59) 

X + 1 


= Nqq{l-ylh{x)) , (63) 

which, in contrast to the bare second law, leads to a 
bound on power. Specifically, this bound is found by 
bounding Lqq in (l5^ using (l63t and then maximizing the 
resulting function with respect to y [s^. This procedure 
yields the simple result 


P(ry,x) < P{r],x) = 4Po 


f?y(l-ry) 

[fj-fj^/x^ 


for 

for 


|x| > 1 

II -I ’ 

X < 1 


(64) 


where we define the standard power 

PO = T.P^qNqqlA (65) 

and y becomes 


y*{x,r]) 

y*{x,-q) 


Axrj 

(x - 1)(1 -I- X - 2?y) 
47y 

x - x^ -2fi + 2xfj 


for |x| > 1, 
for |x| < 1 (66) 


This result is a priori surprising, since this analysis has 
fully incorporated the constraints imposed by the second 
law. 


B. A New Constraint 


We will now prove the existence of an additional con¬ 
straint on the kinetic coefficients (IT^ . which, so far, has 
been missing in our considerations. To this end, we in¬ 
troduce the symmetric matrix 


(Nq, 

. = I Lfj'i 


L 


m 


{L 


wq "t Pqw') 


2 (,Pwq "t Pqw') 


L 


m 


(60) 


through the optimization. Remarkably, for any |x| > 1 
the bound (1641) . which is our third main result, restores 
the quadratic relation between power and efficiency (1581) , 
which we found in our first analysis not invoking the new 
bound dMl) only for the symmetric case x = 1. Con¬ 
sequently, we have shown that, in the linear response 
regime, the power of any cyclic heat engine comprised by 
our theoretical framework must vanish at least linearly as 
its efficiency approaches the Carnot value. We emphasize 
that this quite natural result can neither be derived from 
the laws of thermodynamics nor from micro-reversibility, 
which appears in form of the reciprocity relation (1421) . 
Instead, it relies on the additional constraint (15^ , which 
is beyond both of these principles. 


where 


IV. AN ILLUSTRATIVE EXAMPLE 


Nqq^~i5gqL^^5gq)) 


(61) 


A. Model and Kinetic Coefficients 


will play the role of a normalization constant and ((•)) was 
defined in (1251) . The matrix A has the nontrivial property 
of being positive semidefinite such that the determinant 
of any of its principal submatrices must be nonnegative, 
as we show in appendix[C]by using only the rather general 
assumptions of Sec. ttm 

Two important implications follow immediately from 
this insight. First, by taking the determinant of the lower 


A particularly simple setup for a stochastic heat engine 
consists of a Brownian particle in one spatial dimension 
confined in a harmonic potential of variable strength /c(t) 
and immersed in a heat bath of time-dependent temper¬ 
ature T{t) as schematically shown in Fig. [3] Originally 
proposed in [l3 |. this model has recently been realized 
in a remarkable experiment 0 and can be used to illus¬ 
trate various aspects of stochastic thermodynamics like 
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FIG. 3. Operation cycle of a Brownian heat engine. The 
vertical axis corresponds to the normalized time-dependent 
strength of the harmonic trap in units of p = K(f)/Ko, the 
horizontal axis to the normalized width V = {x^)I(2xq) of 
the distribution function. This plot is analogous to the pres¬ 
sure volume diagram of a macroscopic heat engine such that 
the area encircled by the colored lines quantifies the work 
extracted per operation cycle. Specifically, the plots were ob¬ 
tained using the protocols dZSl) and (| 80 ll for different values 
of the shape parameter d, 2 ^koT = 1, rjc = TcJ-q = 1/10 and 
77 = 1/2. The small graphics show sketches of the potential 
(gray line) and the the phase space distribution, whose color 
reflects the temperature of the heat bath, at the respective 
edges of the cycle. For further explanations, see Sec. nYi 


the role of feedback [H, |4^ and shortcuts to adiabaticity 

M- 

Here, by applying our general theory developed in the 
last sections we will calculate the kinetic coefficients for 
this stochastic heat engine and optimize the protocol K{t) 
controlling the trap strength to obtain maximum power 
for given efficiency. 

In the overdamped limit, due to the absence of kinetic 
energy, the Hamiltonian of the system 

H{x,t) = Ho{x) + AHgyj{x,t) with 

2 

Ho(x) = y^^^ AH=kxI, g^{x,t) =(67) 

depends only on the position x of the particle. Here, kq 
is the equilibrium strength of the trap, k quantifies the 
strength of the time-dependent driving, ^w{t) denotes 
the driving protocol and xo = \/2kBTcl kq is the charac¬ 
teristic length scale of the system. The time evolution of 
the probability density p(x, t) for finding the particle at 
the position x is generated by the Fokker-Planck operator 

L(t) = g.{Ko + Kgnj{t)) dq:X + iikBT(t)dl, (68) 



FIG. 4. Plots of the temperature protocol 'yq{t,d) defined in 
(ISOp (dashed lines) and the corresponding optimal protocol 
'^^{t,g,d) for the trap strength (solid lines) obtained from 
([Till for four different values of the parameter d as functions 
of t/T [13]. For all plots, we have set 2fiK.oT = 1 and fj = 1 / 2 . 
Additionally, the optimal protocol 7 m(t) has been rescaled by 
the dimensionless factor Kogclf^- 


which, in equilibrium reduces to 

L° E gKodj;X + nkBTcdl- (69) 

Here, denotes the mobility and the temperature T(t) = 
T’cT’h/(Th - ^T^q{t)) oscillates between the cold and the 
hot levels Tc and Th = Tc + AT. The equilibrium fluctu¬ 
ations showing up in (jiH)) read 

5ga{x,t) = ^a{t)Koia - kBTjKo) (70) 

with 


= l/(4fcBTc) and ^q = -l/2. (71) 

Formula (liUI) can be easily evaluated by using the de¬ 
tailed balance relation (1321) to transform L° into since 
it is readily seen that the function x^ - kBT^I kq is a right 
eigenvector of with corresponding eigenvalue -2/rKo- 
The resulting kinetic coefficients 


LolP ~ 




T 




- (72) 


are functionals of the protocols 7a (^)- Note that, be¬ 
sides the general reciprocity relation (H^ . these coeffi¬ 
cients also satisfy the special symmetry relation (03]), 
since the factorization condition (1441) is fulfilled in the 
example discussed here. 


B. Optimization 

The optimal protocol 7(),(t,??) for the strength of the 
harmonic trap for a given time dependence of tempera¬ 
ture 7 q(t) maximizes the power output at fixed efficiency 
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rj. It is determined by the variational condition 


0 = 


hw{t) 




(73) 


PIJ^=r} 


where the power P and the heat flux Jq are regarded as 
functionals of 7 u,(i) and 7 q(t)- As we show in appendix 
El this constrained optimization problem has the general 
solution 


A particular advantage of our approach is that it al¬ 
lows to treat situations with a continuously varying tem¬ 
perature of the environment on equal footing with the 
scenario proposed in 1 ^, which involves instantaneous 
switchings between a hot and a cold reservoir. In order 
to illustrate this feature, we consider the specific choice 


, ,, Vl + fisin[27rt/T] 1 

7„ (t, d) E — , -I- - 

2'\/sin^[27rt/T] + d 2 


(80) 


- Vhq{t) 

- 2fjfJ,Ko ^dr(7,(T)-7q))-1-70 (74) 

for 0 < t < T, where we used the abbreviations (1551) , 

%=^f^dtjq{t) (75) 

and 7 o denotes an arbitrary constant. Using this proto¬ 
col, the maximum power 

-Pmax(??) = - J^dt {'^q{t)-Jq) , (76) 

can be extracted per operation cycle at efficiency 77 . 

C. Comparison with General Bound 

In order to compare this result with the general bound 
(IMl) , we evaluate the normalization constant 


for the protocol 75 (t), which, in the linear response 
regime, is proportional to the temperature T(t). The 
function (l80l) . which interpolates between a step func¬ 
tion (d ^ 0 ) and a simple sine (d ^ 00 ) [^, is plot¬ 
ted together with the corresponding optimal protocol 
7 j),(t,r 7 ,d) [5l| in Fig. |4l We find that, for d = 0, this 
protocol shows two sudden jumps occurring simultane¬ 
ously with the instantaneous changes of the bath temper¬ 
ature. Such discontinuities were shown to be typical for 
thermodynamically optimized finite-time protocols con¬ 
necting two equilibrium states [H^. Since, here, we are 
concerned with periodic states generated by permanent 
driving rather than a transient process with equilibrium 
boundary conditions, it is, however, not surprising that 
both, the temperature and the trap strength protocol, 
become continuous as the shape parameter d is increased 
and, respectively, converge to two sines in phase in the 
limit d ->■ 00 . 

For the protocols (1751) and (150]) . the parameter ip be¬ 
comes 


2(1 -r d) -t- ■^d( 1 -t- d) 


(77) 

defined in (I5T|) and rewrite ((75|) as 

-Pmax(?7) = ‘i'^pPovil - fj), (78) 

where Pq is the standard power introduced in (1641) . The 
dimensionless factor 


0 <ip = 


J^dt {lq{t)-lqy 


ifitlqit) 


< 1 


(79) 


quantifies how close the maximum power found in the 
optimization comes to the general bound (IMl) . Since, 
0 < 7 g(t) < 1, it is reached only for 7 , ^ 0. This limit, 
however, requires 75 (t) ->• 0 for any t and thus, inevitably, 
leads to vanishing absolute power. 

For an illustration of this issue, we chose 7 q(t) as the 
step function (1571) such that the system is alternately in 
contact with a hot and a cold bath, respectively during 
the time intervals 7i and T - 7i. We then find '0=1- 
Ti/T and Pq -- 7i/T. Thus, as 71 is decreased, ip comes 
arbitrarily close to 1 and Pq decays linearly to zero. This 
example shows that our bound is asymptotically tight. 


This function decays monotonically from 1/2 for d = 0 to 
1/3 for d ->• 00 . Consequently, ip can not reach its maxi¬ 
mum 1 within the class of protocols (l80l) . This limitation 
can be understood from the argument given below dzg), 
since, for any d, we have 7 ^ = 1 / 2 . 

The standard power Pq is proportional to the function 
2 H-d- + d ), which decays monotonically from 2 to 

3/2 as d increases from 0 to infinity. Thus, Pq exhibits 
the same qualitative dependence on the shape parameter 
d as the efficiency. We can therefore conclude that, at 
least within the model considered here, a steeply rising 
and falling temperature performs better than a smoothly 
changing one. 


V. CONCLUDING PERSPECTIVES 

In this work, we have demonstrated that non¬ 
equilibrium periodic states, which emerge naturally in 
periodically driven systems, can be endowed with the uni¬ 
versal structure of irreversible thermodynamics. More¬ 
over, by using a quite general stochastic approach, we 
have proven the generalized reciprocity relation (1421) for 
the kinetic coefficients characterizing the linear response 
regime. 















10 


Our new framework is particularly useful for a sys¬ 
tematic study of the performance of cyclic heat engines. 
Within the linear regime, bounding the power of these 
machines in such a way that the rather peculiar option 
of Carnot efficiency at finite power is ruled out requires 
the new relation (l62|l . This constraint is beyond the laws 
of thermodynamics and the principle of microscopic re¬ 
versibility and has been proven here for the first time 
on a general level. Remarkably, up to the normaliza¬ 
tion factor, an identical bound has been discovered only 
recently in a numerical analysis of a particular class of 
thermoelectric heat engines [39| . Whether this intriguing 
similarity suggests the existence of a so far undiscovered 
universal principle that applies to periodic as well as to 
steady states leading to a bound on power for any type 
of heat engine that operates in the linear regime remains 
an exciting topic for future research. 

A promising starting point for investigations in this 
directions might be found in the Green-Kubo relations, 
which follow from first principles and provide general ex¬ 
pressions for the conventional kinetic coefficients in terms 
of equilibrium correlation functions [d^. Using a Fokker- 
Planck approach, we have shown that an analogous repre¬ 
sentation for the kinetic coefficients exists in periodically 
driven systems. The quantities related by the relevant 
correlation functions are, however, well defined irrespec¬ 
tive of specific dynamics governing the time evolution of 
the phase space density. It might therefore be possible to 
obtain Hamiltonian-based expressions also for the peri¬ 
odic kinetic coefficients introduced in this work. Finding 
a proper way to take the time dependence of temperature 

I 


into account is arguably the major challenge here. 

This problem also prevents an immediate extension of 
our formalism to the quantum realm. While the first 
part of our analysis, the identification of proper fluxes 
and affinities, carries over to quantum mechanics line by 
line, it is not clear at the moment whether the constraints 
on the kinetic coefficients obtained here classically can be 
likewise transferred or properly generalized. This topic 
appears all the more urgent in the light of recent de¬ 
velopments iiHil showing that the emerging field 
of quantum thermodynamics nowadays comes within the 
range of experiments. 

Appendix A: Derivation of the Expression (1 40 II for 
the Onsager Coefficients 

For an expression of the kinetic coefficients (ITOl) de¬ 
pending only on equilibrium quantities the perturbations 
L^(t) showing up in the linear response solution (IMl) 
have to be eliminated. To this end, we invoke the prop¬ 
erty (1571) of the full Fokker-Planck operator. Substituting 
(1551) into (1571) . expanding the exponential in AH and AT 
and collecting linear order terms provides us with the re¬ 
lations 

where we used the definition (HU). Up to corrections of 
order A^, the periodic distribution can thus be 

rewritten as 


7- r 0 

/(x,t) =p®‘i(x)-I- XI / dr s'-^L°ga{x,t-T)p'^^(x) (A2) 

a=w,q -'0 

7- r 0 

= p^^iK)+ Y, ^ / dTe^^L%^(x,t-T)p^^(x) (A3) 

a=w,q -70 

= p'^^{x)- Y ^(sga{^,t)p'^^{x)- fliTe^°^6Ux,t-T)p<^^x)). (A4) 

a=w ,q V do J 


In the second line, we replaced ga{x,t) with its equilib¬ 
rium fluctuation defined in (1261) . This modification does 
not alter the right hand side of (|A2I) . since L°p®‘^(x) = 0. 
However, it ensures that the function, on which the ex¬ 
ponential operator acts in the third line, which was ob¬ 
tained by an integration by parts with respect to t, has 
no overlap with the nullspace of L° and thus the integral 
with infinite upper bound is well defined. Note that the 
upper boundary term vanishes, since the operator is 


nonpositive [dlj. Inserting (|A2I) into m and (1T51) yields 
I 1 r°° 

Lap = UaSgp)) dr (( 5 «( 0 );%(-t))) . 

(A5) 

Herein, obviously, ga(x, t) can be replaced by its equi¬ 
librium fluctuation Sga{x,t) in the first term. Since any 
constant lies in the left nullspace of L° and, hence, is or¬ 
thogonal to the function e'“ '^Sga{x,t-T)p^^{x), the same 
replacement can be carried out in the second term with¬ 
out leading to additional contributions such that (I40p is, 
finally, obtained. 
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Appendix B: Reciprocity Relations 


In this appendix, we establish the reciprocity relations 
(gH) and (H^ . To this end, we first recall formula dini), 
which becomes 


1 


d"x |(55a(x,t)(55/3(x,t)p®"‘(x) 
dr dga(x,t)e‘- ’■(5g/3(x,t-r)p®‘i(x)j, (Bl) 


using the definitions (l25|) and dm)- By applying the de¬ 
tailed balance relation (15^ and changing the integration 
variable t according to t ^ T - 1 , this expression can be 
transformed to 


L„; 3 [id(x,t),T(t),B] 
1 


(x,-t) % (x,-t(x) 


J dtj d"x 

dr (5g/3(x,-t)e^ ^(55a(x,T-t)p®‘’(x)j, (B2) 


where we introduced the shorthand notation L° e L°(£x). 
Furthermore, to transfer the variable r from the argu¬ 
ment of 6 gi 3 in m to the argument of Sga in (jB2l) . we 
used the identity 


J dt a{t)b{t - t) = J dta{t + T)b{t) 

= J dt a{t + T)b{t) + J dta{t + T)b{t) 

- f dt a{t + T)b{t) = f dt a(t+T)b(t), 
Jt-t Jo 


(B3) 


which holds for any two T-periodic functions a{t) and 
b{t). Finally, we reverse the magnetic field as well as 
the driving protocols and apply the change of integration 
variables x ^ ex, whose Jacobian is 1. By exploiting 
the symmetry dga{s:x.,t) = 6 ga{:>i.,t), which follows from 
condition (|31l) . and carrying out one integration by parts 
with respect to t in the first summand, we obtain 


= dt jd^y. |d5;3(x,t)dg„(x,t)p®‘^(x) 

+ dr 5gp{yi,t)e'-5ga{yi,t- (B4) 
= L/ 3 a[id(x,t),r(<),B] (B5) 

thus completing the proof of the reciprocity relation (1421) . 


We now turn to the special case where the function 
5 u)(x, t) can be separated in the form 


5u,(x,t) = gu,(x)7u,(t). (B6) 

Plugging this expression into (IBID and invoking the def¬ 
inition gq{yd) = -ido(x) yields the expression 


(^) I I B] 

= dt Jd^x |7c(t)7/3(t)dga(x)%(x)p®"'(x) 

+ JJdT'yait)'ji3{t-T)Sga(x)e'- ^(55;3(x)p®'^(x) j, (B7) 

which, by virtue of the detailed balance relation (I32L 
equals 


Lal3[jw{t '), 75(t), B] 

= ^ ^ia{thi3it)Sga{x)dg0{x)p^^{x) 


- JJdT jait)'ji3(t- T)6g/3ix)e'- '^(55a(x)p®'^(x) j. (B8) 


Relation (H51) can now be obtained by following the same 
steps as in the general case, that is, by applying the trans¬ 
formation X ^ ex, reversing the magnetic field and in¬ 
terchanging the arguments 7 u,(t) and jq{t) of Lajs, using 
the symmetry Sga{x) = 6 ga{£x) and performing one in¬ 
tegration by parts in the first term. 


Appendix C: Positivity of A 

The proof that the matrix A defined in (ISOl) is posi¬ 
tive semidefinite consists of two major steps. First, for 
arbitrary numbers ywiVq s Rj we consider the quadratic 
form 


QoiywjUq) - ^ Lapyayf}, (Cl) 

CK,^=w,q 

= ^ f^dr ({G(0);G(-t))) 

Kb -'0 

= ]^fJdtf^^dT (G(0e'”'"G(t-r)) 

where we used the expression go]) for the kinetic coeffi¬ 
cients, defined 


G{x,t)= Y, yadga{x,t) (C2) 

cx=w,q 

and applied the detailed balance relation (15^ to obtain 
the second line. We recall the definitions (1251) for the 
meaning of the angular brackets. The crucial ingredient 
for this first step consists of the identity 
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2 k 


^ ^ - r))e^°'^'Git - r')) + ((e^^'^Git - T)){dr^e^°'-')G{t - r'))| 


- 7 ” ^ oo ^ oo 

Jo 

1 r'T' 


dt J^llT\2l^G{t)e^°^^G{t-T)) + dt J^hr ((e^“^^G(i - T))e^°^^'G(t - r')) 


2fc_BT -Vo -Vo 
1 r'^. 


j^^j-£dt£ dr l^G(t)e^°^'"G{t-T)'^ = Qo(y^,yq). 


(C3) 

(C4) 

(C5) 

(C6) 


Here, we used the relation 

= (h[°U) , (C7) 

which holds for any functions H(x),_B(x) by virtue of 
the detailed balance condition (j32|) . to obtain (IC4I1 from 
(IC3I) . Expression (IC5I) follows by applying an integration 
by parts with respect to r and r', respectively, in the first 
and the second summand of (EH). Finally, the second 
contribution in (IC5I) vanishes after carrying out the t- 
integration, since the function G(x, t) is T-periodic in 
time. 

Next, we note that (IMl) implies 
(A(L°Ul°1')H) (C8) 

= J d"xH(x)(LV‘’(x) +p®'^(x)L°^)H(x) 

= Jd^x [p®‘i(x)]5H(x)(K° + K°^)[p"'i(x)]5H(x), 
where 

K° E [p®'i(x)]-k°[p®‘i(x)]i (C9) 

Since the Hermitian part of this operator is negative 
semidefinite [4l|, it follows that (IC8I) is non-positive for 
any H(x). Hence, we can conclude that the quadratic 
form Qo{yw,yq) is positive semidefinite, since it can be 
written in the form (1C3|) . 

For the second step of the proof, we introduce the 
quadratic form 

Q{yw,yq,z) = Qo{yw,yq) + Qiiyw,yq,z) (CIO) 

with 

Ql i,yw: Vq 7 '^) = qq^ +22 Lq^yci 

ot=w,q 

= ({{Fl^^F)) - 2 {{FG} - 2fJdT ((E(0); G(-r)))) 

= |2 (F(t)L°tF(t)) -4(E(t)G(t)) 


where yw^Vq, z eM. and 


F{x,t) = z6gg{x,t). (C12) 


Note that, in (IC11|) . we used (HOI) and (ICT) as well as the 
detailed balance condition (l32ll . The expression (jCllIl 
can be rewritten as 

Q,iy^,yq,z) = fjdt |(F(t)(L°t + L°t)F(t)) 

+ JJdT (F(t)(L°+ + L°t)(eL”'"G(t-r))) 

+ ^'dr (e^°^^G(t-T))(L°UL°1')E(t))|. (C13) 

This assertion can be proven by expanding (IC13I) , invok¬ 
ing (Ell) as well as the identity 


= zL°^gq{Ti,t) = -zjq(t)L°'^ Ho{x.) 

= -ziq{t)L°'^ Ho{x) = L°’l'E(x,t), (C14) 


which is implied by condition (1361) . and integrating by 
parts, first with respect to r and then with respect to t, 
respectively, in the second and third term showing up in 
(IC13F Finally, putting together (IC3I) and E13I) leads to 
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Q(?/„,y5,z) = ^ dt |(F(i) + ^'dTe'-“^^G(t-T)](L°'l' + L°'l')(F(t) + ^ dr'G(i - t'))| • (C15) 


The average showing up in this expression is of the form 
(IC8I) and thus must be non positive. Consequently, we 
have Q{yw,yq, z) > 0 for any yw,yq,z. Moreover, since 

Q(yw,yq,z) = y* Ay (C16) 

with y E (z, yw,yqY it follows that the matrix A must be 
positive semidefinite and thus the proof is completed. 


Appendix D: Optimal Protocol 


The aim is to determine the optimal protocol 7 ^(t), 
which maximizes the rescaled output power 


P = 




= -{LwwX^ + LwqX) with x = ^w/d^q (DY 


for given time dependence of the bath temperature gq(t) 
and normalized efficiency 


V = 


Lww\ + LwqX 


(D2) 


AqwX ^qq 

This task is captured by the objective functional 

(A ~ 1)X ^ww (A “ ^'nX.^qw A?]Z/qg, (D3) 

where the additional constraint (imi) is taken into ac¬ 
count by introducing the Lagrange multiplier A. By in¬ 
serting dm]) becomes 


'Pbw{t),Xqit),X] = ^ [ dt Y, 'U'ap\ja{t)'y 0 it) 

‘ ^ OL,0=w,q 

-^“dr7„(t)7;3(<-r)e-2''"'>"|. (D4) 


Here, we introduced the coefficients 


luu,w ""-9 ) = _2fcBT2 ( ' 

\ y^qw Uqq / y 


2/(A-1)x^C (A-1)xC«;6 




A^C, 


) 


(D5) 

for notational simplicity. The convolution type structure 
of (ID4I) naturally suggests to solve the variational prob¬ 
lem by a Fourier transformation. We expand 


Xa{t) = Y 

neZ 

and thus obtain 


a incot 

c„e 


with 


' E 27r/r (D6) 


= {- 2 y,Ko) Y T. '^apCnC-nZ 


inuj 


iriLu - 2fino 


(D7) 


Since (lD7l) is quadratic in the Fourier coefficients c“, it 
is straightforward to carry out the optimization with re¬ 
spect to c“. Taking into account that the protocols must 
be real and therefore c“„ = c“*, the conditions 


yield 


c 


W 

n 


9c^'P[lwit)yq{t),X] = 0, 

dc^,V[ywit),Xq{t),X] = 0 


( ^wq “*■ • '^wq '^qw \ n 

AU'njiv 2iTILUU'hj'hj / 


(D8) 

(D9) 

(DIO) 


for n 0. Note that (ID7I) does not depend on Cg and 
thus the optimal protocol will be unique only up to a 
trivial offset Cq . To comply with the constraint (ID2I) . 
the Lagrange multiplier A must be chosen such that 


dxV[yUt),Xq{t),X] = 0, (Dll) 


where the derivative has to be taken before the cJC are re¬ 
placed by the solution (IDlOp . After some algebra, (|D11I) 
reduces to the simple condition 

f-(A-1)2(7-1)2 = 0, (D12) 

which is fulfilled for 


A± = (l-7±7)/(l-7). (D13) 

Inserting (ID13P and (IDIOI) into (ID7I1 yields 
P^^V[yUt)nqit),X^] = 0, (D14) 

oo 

P- = 'P[lwit),yq{t),X-] = SksT^- fj) Y K\^ 

n=l 

= 17 ( 1 - 1 ?)^ dt iXqit)-%f , (D15) 

where 7 ^ is defined in (1751) . Consequently, the relevant so¬ 
lution for the Lagrange multiplier is given by A_. Finally, 
the optimal protocol x*(t,r]) is obtained by summing up 
the Fourier series (ID 6 I) . The explicit result (l7^ can be 
found by, first, evaluating 


7;(t,7)=*a;^nc-e“-‘ 

neZ 

iTlUj(^U^q + Uqw') 
2 ?/ 

iq{t) + 2AiK0 


-e(- 

neZ V 


- 2/iKo 


■^wq ^qw 


2ll 

‘2k bT^ 


2 ut, 


Uwq Uqw \ ^q ^inujt 

9 I 

{lq{t)-iq) 


X 


(2^K07 {lq{t) - Iq) - (1 - ’7)79(^)) (D16) 


and, second, solving the simple differential equation 
(IDIB . 
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